EXISTENCE OF WEAK SOLUTIONS FOR A NON-CLASSICAL 
SHARP INTERFACE MODEL FOR A TWO-PHASE FLOW OF 
VISCOUS, INCOMPRESSIBLE FLUIDS 



HELMUT ABELS AND MATTHIAS ROGER 

Abstract. We introduce a new sharp interface model for the flow of two immis- 
cible, viscous, incompressible fluids. In contrast to classical models for two-phase 
flows we prescribe an evolution law for the interfaces that takes diffusional effects 
into account. This leads to a coupled system of Navier-Stokes and Mullins-Sekerka 
type parts that coincides with the asymptotic limit of a diffuse interface model. We 
prove the long-time existence of weak solutions, which is an open problem for the 
classical two-phase model. We show that the phase interfaces have in almost all 
points a generalized mean curvature. 



1. Introduction 

We study the flow of two incompressible, viscous and immiscible fluids like oil and 
water inside a bounded domain f2 or in f2 = T d , d = 2, 3. The fluids fill domains £l + (t) 
and Q~(t), t > 0, respectively, with a common interface T(t) between both fluids. 
The flow is described in terms of the velocity v : Q x (0, oo) — > IR d and the pressure 
p: Q x (0, oo) — > R in both fluids in Eulerian coordinates. We assume the fluids to be 
of Newtonian type, i.e., the stress tensors are of the form ^(v^p) = v ± Dv — pI with 
constant viscosities > and 2Dv = Vi> + Vv T . Moreover, we consider the case 
with surface tension at the interface and assume that the densities are the same (and 
set to 1 for simplicity). For the evolution of the phases we take diffusional effects into 
account and consider a contribution to the flux that is proportional to the negative 
gradient of the chemical potential ji. Precise assumptions are made below. 

To formulate our model we introduce some notation first. Denote by n the unit 
normal of F(t) that points inside Q + (t) and by V and H the normal velocity and scalar 
mean curvature of T(t) with respect to n. By [•] we denote the jump of a quantity 
across the interface in direction of n, i.e., [f](x) = \im h ^ (f (x + hn) — f(x — hn)) for 
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x G r(i). Then our model is described by the following equations 



— divT ± (f ,p) 


= 


in ^(t),* > 0, 


(i.r 


divv 


= 


in tt ± (t),t> 0, 


(1.2) 




= 


in fi^t),* > 0, 


(1.3) 


-n ■ [T(v,p)] 


= kHh 


on r(t),t > 0, 


(1.4) 


V 


= n ■ v\r(t) ~ m l n " V/i] 


on r(t),t > 0, 


(1.5) 


^|r(t) 


= kH 


on r(t),i > 0, 


(1.6) 


v\an 


= 


on dQ, t > 0, 


(1.7) 


n n ■ mV//|an 


= 


on <9fi, t > 0, 


(1.8) 


n + (0) 


= ^o + 




(1.9) 


«|t=Q 


= vo 


in f2, 


(1.10) 



where Vq, are given initial data satisfying 8Qq D dQ = and where k, m > are 
a surface tension and a mobility constant, respectively. Implicitly it is assumed that 
v,fi do not jump across T(t), i.e., 

[ v ] = [ fl ] = onr(t),t>0. 



Equations ( 1.1 )-( 1.2 ) describe the conservation of linear momentum and mass in both 



fluids and (1.4) is the balance of forces at the boundary. The equations for v are 



complemented by the non-slip condition (1.7) at the boundary of fi. The conditions 



(1.3), (1.8) describe together with (1.5) a continuity equation for the mass of the 
phases, and (1.6) relates the chemical potential [i to the L 2 -gradient of the surface 
area, which is given by the mean curvature of the interface. In this formulation of 
the model we assume (if Q ^ T d ) that F(t) does not touch dQ. 

For m = the velocity field v is independent of /x. In this case, (1.5) describes 



the usual kinematic condition that the interface is transported by the flow of the 



surrounding fluids and ( 1.1 )-( 1.10) reduces to the classical model of a two-phase 
Navier-Stokes flow as for example studied by cf. Denisova and Solonnikov j6], where 
short time existence of strong solutions is shown. On the other hand, if m > 0, the 
equations (1.3), (1.6), (1.8) with v — define the Mullins-Sekerka flow of a family of 
interfaces. This evolution describes the gradient flow for the surface area functional 



with respect to the H l (Q) scalar product. Therefore we will also call (1.1 )-( 1.10 ) 
the Navier-Stokes/Mullins-Sekerka system. 



The motivation to consider ( 1.1 )-( 1.10) with m > is twofold: First of all, the 
modified system gives a regularization of the classical model m = since we change 
from a parabolic-hyperbolic system to a purely parabolic system (cf. also the effect 
of m > in ( 1.13[ ) below). Secondly, (1.1)-(1.10) appears as sharp interface limit of 
the following diffuse interface model, introduced by Hohenberg and Halperin [9] and 
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rigorously derived by Gurtin et al. [8]: 



div(u(c)Dv) + Vp = 


-e div(Vc ® Vc) 


in x 


(0,oo), 


(1.11) 


divv = 





in x 


(0,oo), 


(1.12) 


d t c + v ■ Vc = 


mA/i 


in x 


(0,oo), 


(1.13) 


^ = 


e" 1 /' (c) - £ Ac 


in x 


(0,oo), 


(1.14) 


v\an = 





on <9f2 


x (0,oo), 


(1.15) 


d n c\aQ. — d n ^\an = 





on <9f2 


x (0,oo), 


(1.16) 


(v,c)\ t=0 = 


(«o, c ) 


in Q. 




(1.17) 



Here c is the concentration of one of the fluids, where we note that a partial mixing 
of both fluids is assumed in the model, and / is a suitable "double-well potential" 
e.g. f(c) = c 2 (l — c) 2 . Moreover, e > is a small parameter related to the interface 
thickness, \i is the so-called chemical potential and m > the mobility. We refer to 
PQ for a further discussion of this model and to the appendix where we discuss the 
convergence of ( 1.11 )-( 1.17 ) to varifold solutions of (1.1 )-( 1.10 ). 



Sufficiently smooth solutions of ( 1. 1 )-( 1.10 ) satisfy the following energy equality 



\T(t)) = - / u(X)\Dv\ 2 dx - m / \Vn\ 2 dx, (1.18) 



where 2/(0) = V- and u(l) = u + and 7i d ~ 1 denotes the (d — l)-dimensional Hausdorff 
measure. This identity can be verified by multiplying ( 1.1[ ) and (1.3) with v, ft, resp., 
integrating and using the boundary and interface conditions (1.4)-(1.8). This energy 



equality motivates the choice of solution spaces in our weak formulation and shows 
that the regularization introduced for m > yields an additional dissipation term. 
In particular, we expect fi(-,t) G H 1,2 (Q) for almost all t G M+ and formally, using 
Sobolev inequality and ( 1.7[ ), that H(-,t) G L 4 (r(t)) for d < 3. This gives some 
indication of extra regularity properties of the phase interfaces in the model with 
m > 0. 



Our main result is the existence of weak solutions of (l.l)-(l.lO) for large times. 
For the definitions of the function spaces we refer to Section 2] below; the concept of 
generalized mean curvature for non-smooth phase interfaces is taken from [16J, see 
Definition 14.41 below. 



Theorem 1.1. Let d = 2,3, T > 0, let Q C M. d be a bounded domain with smooth 
boundary or let Q = T d and set Qj< = Q x (0,T). Moreover, let z/(0) := v_, 
i/(l) := z/ + and K,rn > 0. Then for any v G L 2 (fi), X G BV(fl; {0, 1}) 
there are v G L°°(0, T; L 2 (O)) n L 2 (0, T; H^(Q) d ), X G L°°(0, T; BV(Q; {0, 1})), 
// G £ 2 (0, T; i? 1 (0)) ; i/iat satisfy (1.1)-(1.10) in the following sense: For almost 
all t G (0, T) the phase interface d*{X(-,t) = 1} has a generalized mean curvature 
vector H(t) G L s (d\V X h {t)\) d with s = 4 if d = 3 and 1 < s < oo arbitrary if d = 2, 
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such that 



/ (-vd t <p + v ■ Vv<p + u(X)Dv : Dip) d(x, t) 



ip\t=o-v dx= k / H(t)-ip(t)d\VX(t)\dt, 
< > Jo Jn 



1.19) 



holds for all <p G C7°°([0, T}; C%> (fi)) wt/i y? 



#(a^ + div(^i;))da:d* + / A'o(x)^(0,x)rfx =m / V^-V^dxdt (1.20) 



fco/ds /or a/I ij) G l7°°([0,T] x Q) with ip\ t=T = and 
V*(-,t) 



^•) = A*(*,v, v ^| ( . jt) 
ZioZds /or almost all < t < T . 



H — almost everywhere on d*{X(t, .) = 1} (1.21) 



Remark 1.2. (1.19) is the weak formulation of (1.1), (1.4), and (1.10). It is obtained 
from testing (1.1) with tp in fi ± (t), integrating over fi + (t) U Vt~{t) and using (1.5) 



together with Gauss' theorem. Similarly, (1.20) is a weak formulation of (1.3), (1.5), 
(1.8), and (1.9). The conditions (1.2), (1.7) and [v] — \p] = on T(t) are included 



in the choice of the function spaces, namely v(t) G Hq(Q), n(t) G if x (f2) for almost 
every t > 0, and (1.6) is formulated in ( 1.21[). 

Finally, we note that, because of (1.21), (1.19) is equivalent to 



/ (—vd t ip + v ■ Vvip + v{X)Dv : Dip) d(x, t) 
Jn T 



— / VQ-p\t=odx= — \ XVfi ■ <pd(x,t) 

'U Jq t 



1.22) 



for all ip G C°°([0,r]; C^ a {Q)) with p\ t=T = 0. The latter form will be used for the 
construction of weak solutions. 

Remark 1.3. Compared to the available long-time existence results for the classical 
model m = and as a consequence of the diffusive effects that are included in our 
model, Theorem 11 yields the long-time existence of more regular solutions. In 
the classical case m = Plotnikov [15] and Abels [21 E] have shown the long-time 
existence of generalized solutions. However, in their formulations the phase interfaces 



are in general not regular enough to define a mean curvature. Condition (1.4) is only 



satisfied by a varifold that may depend on the construction process of the solutions, 
that is in general not rectifiable, and lacks a (d — l)-dimensional character (due 
to concentration and oscillation effects of the interface and e.g. the formation of 
"infinitesimal small droplets", cf. the discussion in |3|). In contrast, in our weak 
formulation the phase interfaces have a generalized mean curvature that enjoys the 



integrability property that are expected, in the smooth case, from (1.7), the energy 



equality, and the Sobolev inequality for the chemical potential. A similar result for 
the case m = is an open problem. 
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We note that a similar but different regularization was proposed by Liu and 



Shen [10]. In their model (1.5) is replaced by 

V = n ■ v\r(t) + mH. 

Local in time well-posedness for the latter system was proved by Maekawa [T2] 
Physically, this model has the disadvantage that the mass of the fluids, i.e., (^(t)) 



is not preserved in time, while this is the case for our system (1.1 )-( 1.10 ). 



Remark 1.4. We note that our concept of weak solution does not include a formu- 
lation of a contact angle condition in the case that Q is a bounded domain and the 
phase boundary d*{X(t, .) = 1} meets the boundary of the domain dfl. Even for 
weak solutions of the Mullins-Sekerka flow as constructed in [16] the formulation of 
boundary conditions is an open problem. 

For simplicity we will assume k = m = 1 in the following. All statement below 



and the proof of Theorem 1.1 will be valid for general m, k > if modified accord- 
ingly. The structure of the article is as follows: First the basic notation and some 
preliminaries are summarized in Section [2] Then weak solutions of a time-discrete 
approximate system are constructed in Section |3j Our main theorem is proved in 
Section [4] by passing to the limit in the approximate system. Finally, in the appendix 



we prove the convergence of the diffuse interface model (1.11 )-( 1.17) to (l.l)-(l.lO) 



However, in this limit we have to work with a weaker notion of generalized solutions, 



compared to the notion of solutions that we use in Theorem |1.1 

2. Notation and Preliminaries 



For A,B £ R dxd we denote A : B = tr(AB) and \A\ = VA : A. Given a G l d 
we define a <S> a £ R dxd as the matrix with the entries ajOy, i,j = l,...,d. The 
space of all fc-dimensional unoriented linear subspaces of R d is denoted by G k - If X 
is a Banach space, X* denotes its dual and (x*,x) = (x*,x)x*,x, x* £ X*,x £ X, 
the duality product. If if is a Hilbert space, then (-, •)# denotes its inner product. 
Moreover, we use the abbreviation (■, -)m = (■> -)l 2 (m)- 

For s > we denote by [s] the integer part of s and for / : R — > X we define the 
backward and forward difference quotients by 

_ f(t)-f(t-h) + _ f(t + h)- f{t) 
<>t,hJ - h > a t,hJ - h 

2.1. Measures and BV-Functions. Let X be a locally compact separable metric 
space and let Cq(X; M. m ) by the closure of compactly supported continuous functions 
/: X — > lR m , m £ N, in the supremum norm. Moreover, denote by Ai(X;M. m ) 
the space of all finite R m - valued Radon measures, M{X) := M(X;R), and Mi(X) 
denotes the space of all probability measures on X. By the Riesz representation 
theorem M(X; R m ) = C (X; R m )*, cf. e.g. d Theorem 1.54]. Given A £ M(X; R m ) 
we denote by |A| the total variation measure defined by 

{oo oo 
\ \{A k )\ : A k £ B{X) pairwise disjoint, A = (J A k 
k=0 k=0 
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for every A G B(X), where B(X) denotes the a-algebra of Borel sets of X. Moreover, 
A : X — > M. m denotes the Radon- Nikodym derivative of A with respect to |A|. The 
restriction of a measure fj, to a /i-measurable set A is denoted by (p,[A)(B) = n(AnB). 
Finally, the s-dimensional Hausdorff measure on IR d , < s < d, is denoted by 7i s . 
Let U C R d be an open set. Recall that 

BV(U) — {/ 6 ^(C/) : V/ G A<(C/; 

||/||BV(i7) = ll/IU 1 ^) + ||V/||x(a;Rd), 

where V/ denotes the distributional derivative. Moreover, BV(U; {0, 1}) denotes 
the set of all <Y G BV(U) such that A'(x) G {0, 1} for almost all x Ell. 

A set E C [/ is said to have finite perimeter in £7 if ^ G BV(U). By the structure 
theorem of sets of finite perimeter |VA^| = 7i d ~ 1 \_d*E, where d*E is the so-called 
reduced boundary of E and for all cp G Cq(U, Mr) 

-{VX E ,tp)= I div^dx = - / Lp-n E dU d ^, 

JE Jd*E 

where n E (x) = ly^ji , cf. e.g. Note that, if E is a domain with (^-boundary, 
then d*E = dE and n E coincides with the interior unit normal. 

2.2. Function Spaces. As usual the space of smooth and compactly supported 
functions in an open set U is denoted by Cq°(U). Moreover, C°°(U) denotes the set 
of all smooth functions with continuous derivatives on U. If X is a Banach-space, 
the X- valued variants are denoted by C$°(U;X) and l7°°(17;X). For < T < oo, 
we denote by Lf oc ([0, T); X), 1 < p < oo, the space of all strongly measurable 
/ : (0, T) -> X such that / G L p {0, T"; X) for all < V < T. 
Furthermore, C^(fi) = G Cq 00 ^ : div<^ = 0} and 

cf. e.g. [20] for other characterizations and properties of L^(fi). 
Finally, we will use the following notation: 

Hf 0) (n) : = H\n)n{fi-. [ » = o}, 

W (mo) (Q;{0,l}) := W(fi;{0,l})n{jf Ar = m }, 

fl^(n) := ff? 0) (n)*. 

Here we note that H^(Q) is equipped with the norm H/Hija (n) = ||V/||i,2(n) and 
iif^(fi) with the dual norm associated to the latter norm. In particular, this yields 
the useful relation 

Wf\\ H -Un) = HVC-A^-VH^n) for all / G HzHn), (2.1) 
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where — : HhJfi,) — > HLJQ)* is the weak Laplace operator with Neumann 
boundary conditions defined by 

(-A N w,(p) H i o){n) * tH i o){n) = / Vw-V^dx for all cp G H^Q). (2.2) 

3. TIME-DISCRETE APPROXIMATION 

In this section we will construct weak solutions of an approximate time-discrete 
system. Fortunately the coupling of the Navier-Stokes to the Mullins-Sekerka system 
can be treated explicitely. The main result of this section is: 



Proposition 3.1. Let the assumptions of Theorem 1A_ be valid and let mo = 
j Q Xq dx. Then for all h > sufficiently small there exist time-discrete solutions v h G 
L°°(-h,T;H^(ty), X h G L°°(—h,T; BV( mo )(Q; {0, 1})), fj* G £°°(0, oo; #£ 0) (fi)), 
and Lagrange multipliers X h : [0, T] — * K such that for all t G (0, T) the following 
equations hold 

1 (v h (t) - v h (t -h))r) + v h (t - h) ■ Vv h (t)ri) dx 



+ / u{X h (t))Dv h {t) : Drjdx= - [ X h (t)V^{t) ■ r}dx (3.1) 
Jn Jn 

for allrje H^Q), 

v h (t) = v Q , X h (t) = X (3.2) 

for -h<t<0, 

J ( Xk{t) " f (t ~ k h ~ (v h X h )(t - h) • ve) dx= - JvrtW-VZdz (3.3) 
for all £ G # x (fi), and with /i h (t) := $(t) + X h (t) 

for all 7] G C l (Vt\ IR d ) with i] ■ uq = on <9f2. Moreover, we have the estimates 

sup h\At)\\b<p) + sup / d|V* fc (f)| 
te(o,T) ^ te(o,T) 

+^min||^||i 2( n T ) + \\\V4\\h(n T ) < C(\\V \\ L 2 {S1) ,J^ d|V# |) (3.5) 

and for all t G (0, T) 

IIa*o(*) + A ft (t)||^) < C(l + \\V4(t)\\ LHn) ) (3.6) 

holds, where C depends only on d, Q, T, m , the initial data, and v min = 
min(z/(0),i/(l)). 



In the remainder of this section we prove Proposition 3^ The first step gives the 
solvability and estimates for a time-discrete and regularized Navier-Stokes equation. 
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Lemma 3.2. Let v G H^^n), X G BV {mo) (Q; {0, 1}), fi G H 1 ^), and h > 6e 
given. Then there exists a solution v G H^lfl) of 

' 'v-v)+v- Vv J <pdx+ [ u(X)Dv : Difdx = - f XV 'fi ■ <p (3.7) 



/or a// G Co^(fl). For eac/i solution v G ^^(n) 

^lkll! 2 (fi) + ^ / v(X)\Dv\ 2 < h\v\\h {n) -h I XVfi-vdx (3.8) 

«/ fi */ fi 



holds. 



Proof. We show the existence of a solution t> G Hoo-(^) °f (3-7) that satisfies (3.8) 
with the aid of the Leray-Schauder principle, cf. e.g. [201 Chapter II, Lemma 3.1.1]. 
To this end, we define L: ^Jfi) -> and G: -> £ L^ /2 (ft) 

by 



fi 



y ^— X V/i — — (u — v) — v ■ Vim ■ if) dx 



for all v,<p G ^^(Q), ^ G L^(fi). By the Lemma of Lax-Milgram, L: (SI) — >■ 
ffQ ff (Q)' is an isomorphism. Moreover, it is easy to check that G: HQ a (Q) — > L^(fi)* 
is a continuous mapping, where we note that 5 • Vi> G L 3 / 2 (fi). Since L^(SV)* > 
ifoo-(^)* compactly, G: H(j a (Sl) — > 1/^(0)* is completely continuous. Thus F 



L 1 G f : iJg CT (S7) — ► i?oo-(^) is completely continuous and (3.7) is equivalent to the 
fixed-point problem v = G(v). In order to apply the Ler ay- Schauder principle, let 
R > be such that R 2 = M 2 \\ V/x||| 2(n) + f\\v\\l2 (n) , where M is a constant such 
that 

\M 2 H i(n) <M I v{X)\Dv\ 2 dx for all v G F^fi). 
Jn 

It remains to show that for all v G Hq^Q) 

v = \F(v),\e[0,l] \\v\\ H i A n) < R. (3-9) 

To this end let v = XF(v) for some A G [0, 1]. Then Lv = \G(v) and therefore 

/ (v — v) (p dx + / Xv ■ Vf ip dx + / v(X)Dv : Dip dx = X X\7fi-(pdx 
Jn Jn Jn Jn 



X 
h 



for all ip G Hq^SI). Choosing ip = v yields 

X n i,n f o , . , f , A 



^-|HI^(fi) + / ^(^)|-D^| 2 dx < X I XVfi-vdx+^-WvWl^, (3.10) 

»/ fi */ fi 
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where we have used (v — v) • v 



\v\ 2 — l\v\ 2 + \\v — v\ 2 > hv\ 



using that \J Q XVfi-v\ < || V/i|| L 2 (n) |M| L 2 (n) , 

\M 2 m (n) < M J v{X)Dv: Dvdx<M (\\Vn\\ L 2 



(n)\\v\\L 2 (n) + 



^M 2 . Hence, 



2^ll«lli»(o) 



< — va* j 



2 H • r-ilL 2 (C) + 2^ W' v \\l 2 (q,) ■ 2 

and therefore \\v ||fl-i(n) < R- 

Because of (3.9 ), the Leray-Schauder principle implies the existence of a fixed point 
v G #o CT (fi) which solves dO. Finally, (Q follows from (|3.10b with A = 1. □ 



M 



+ -\\v\ 



Next we solve the appropriate versions of the Mullins-Sekerka part (3.3), (3.4). 
We follow Luckhaus-Sturzenhecker [TT] and use that the Mullins-Sekerka flow is the 
-f/" _1 -gradient flow of the surface-area-functional. 

Lemma 3.3. For X G BV {mo) (Q; {0, 1}) and v G H^(Q) there exist X G 
BV( mo )(fl; {0, 1}), /io G HLJtt), and a constant A G R suc/i £/iai 



V/i • V£ + - (# - X)£ - vX ■ Vn dx 



(3-11) 



/or a// £ G // 1 (fi), snc/i t/iat 



jf (divi7-|||r.I?i7T|^|)d|VAr| = J^X dw ((fi + X) v ) dx (3.12) 



/or a// 77 G C 1 (fi, IR d ) with rj ■ = on <9fi ; and stzc/i that 



d|V*| + -||V/i<,|| 



< 



d|V#| + 



h, 



7~,l| 2 



Moreover, we have 



|A| < C(m ,n,fi)(l + J d\VX\^J d\VX\ + ||V/z 
|/*o + A||fli ( n) < C(m ,n,fi)(l + ^ d|V#|)(jf d\VX\ + \\V/i, 



0||L 2 (Q) y , 

olU 2 (n) 



(3.13) 

(3.14) 
(3.15) 



Proof. We split the proof into several steps. 

Step 1 : There exist /i G flf 0) (fi),Af G BV (wo) (fy {0, 1}) satisfying (gll), gig , 
and enjoying a minimizing property from which we will deduce (3.12). 
Define a functional F h : W (mo) (ft; {0, 1}) — ► R, 



IV 



•^l + wr\W - X + hv ■ VX\\ 2 M _i, , 
2/i" IIH ( )( n ) 



(3.16) 



for a G W( mo )(0; {0, 1}). We remark that v ■ VX G H^(Q)* is defined by 

(iJ.V^,C> = - / Xv-VCdx, (3.17) 
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where we note that divt> = 0. Because of (2.2), 



\a- X + hv- VX\\* 1(m 



[a -X + hv- VX, (-An)- 1 (a -X + hv- VX)) x 



H{ 0) (n)*,Hl 0) (n) 



. (3.1? 



By the L 1 (f2)-compactness of bounded sequences in BV(Q), the lower semi-continuity 
of the perimeter under L 1 (fi)-convergence, and the continuity of the embedding 
L 2 JQ) — > H^^Vl) there exists a minimizer X 6 5V( mo )(fi; {0, 1}) of F h . More- 
over, 



/'o 



-A^^X -X) +v- VX^j 



(3.19) 



satisfies ( [3~TT] ). We deduce now fl3~13| from F /l (A') - F ft (Af) < and (Q, ( |3~l9| . 
In fact, 



Vf-A^r 1 ?! • VX\ 2 dx 



[ d\VX\ + ^- [ \hVfi \ 2 dx < [ dlVXl + ^/v-VX^-A^-^v-VX)) i 

Jf7 Jf7 ^ N ' H (0)' H (0) 

d\VX\ + 
d|VAT| 

where in the last step we have used that by ( 2.2[ ), (3.17) 

|2 

ll, 2 (n) " 



< 



IVI-A^)- 1 ^ V^| 12 



I fill 2 



v-VXA-An^v-Vx) 



/ cj-l fj-1 

•"(O'^to) 



Xv -W(-A N y l v -WX dx 



< ||u||ia(n) II V(— Ajv) £ • V^|| L 2 (n) , 
hence || V(-A A r)" 1 {; • V^|| L 2 (n) < HfJHi^n)- 

Step 2: We compute the first variation of F h in with respect to volume preserv- 
ing variations. With this aim we consider a smooth family (<E> s ) se (_ E)e ) of smooth 
diffeomorphisms $ s : Q — > Q with <3>o — Id and variation field r\ such that 



V 



d_ 
ds 



$ s e C°°(fl,'. 



rj ■ rin = on dQ. 



s=0 



Assume that the variations $ s conserve the volume of {X = 1}, that means that 
cr s := X o satisfy 



a s dx = m for all — e < s < e, 



in particular a s G BV( mo )(Q; {0, 1}) and J n X divrj = 0. Since X minimizes F h in Z 
we have j^\ s= oF h (X s ) = 0. The first part of F h is given by the perimeter-functional 
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Vn and 

d 
ds 



-Me*) = 6V a (X)(T}) = J (divr? - ■ 0»7j|§r) d\VX\. (3.20) 



Let K denote the second part of F . Since (— Ajv) is linear and symmetric, since 

d_\ 
8s \ 



d 

ds 



ocr s = — V X ■ i], and by ( |3.19[ ) we obtain that 
K h (a s ) 



s=0 



5K h {X){rj) = (VX-^^h-i^ = - I Xdiv( Wo )dx. (3.21) 



We therefore deduce from the minimality of X that 







divr/ 



VX 

\vx\ 



Dr) 



VX 

\vx\ 



d\VX\ 



X div(?7/xo) dx. 



(3.22) 



Step 3 : We next prove ( |3~l2| . Fix £ G C °°(^,R d ) such that 

X div£ dx 7^ 0, £ • Tin = on <9fi, 



and choose a family (/i r )re(-ei,ei) of smooth diffeomorphisms of Q with /i = Id 
£. Similarly, for given 77 G C°°(fi,lR d ) with 77 • no = on <9J1 
we let (<7 s ) se (_ £l)£l ) be a family of smooth diffeomorphisms of f2 with go = Id and 



and 



ill h 

Q r |r = 0"T 



d_ 
ds I 



--o9s 



rj. Then the function / : (— e 1 ,£ 1 ) z 



defined by 



f(s,r) 



/ X o (g s o h r ) 1 dx — m$ = / (det ((Dg s o h r )Dh r } — l) dx 
Jn Jn 



for sufficiently small E\ > satisfies 

/(o,o) = o, <9 r | 



=of(s,r) 



Xdiv£ ^ 0. 



Since / is smooth we obtain by the implicit function theorem that there exists < 
e < ei and a smooth function g : (— e, e) — ► (— £i,£i) such that 



f(s,g(s)) 
We therefore deduce that 

hence 



for all s G (—£,£). 



(3.23) 



of{s,g{s)) = j^divr l + g'(0)dw^)Xdx, 



q'(0) = ~(J Xdiv^dx^J 1 J Xdivridx. 



(3.24) 



By (3.23) the family (g s o h r r g )) s e(-e,e) defines a variation of f2 that conserves the 



volume of X. The corresponding variation field is given by 

— I 
ds 



n 



=0 (g s oh r ) = r] + q'(0)£. 
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Step 2 therefore implies that 
= 



divfj — |rT — x ■ Dfj |rT v| ) d| VAf| — / <Y div(r7/i ) fix 



and yields by (3.24) 



^7x wx \ r 

divi] — lr7 v , l • Drj j^ — { ) d\VX\ — / X div(r]fio) dx 

Jn 



VX\ '\VX\ 

vx „vx 



^ |V*| — |VAT| 



^7 



d\VX\ - / Xdiv(Zn )dx 
n 



A / A'divT] 



(3.25) 



with 



A := 



X div £ (ix 



_ ^div£- — -Z^— Jd|V*|- / ^div(£/i )rfx . (3.26) 

J f2 II II Jii 



This proves (3.12). 



Step 4: Finally we derive (3.15) by choosing a particular £ in Step 3. We adapt 



the proof given in j5]. First we choose a Dirac sequence (<p&)&>o with kernel <p G 
C~(£i(0)), < ^ < 1, and set 



1 



Af 5 := ^ * <£} 5 , : 

Let ip : Q — > K be the solution of 

A^> — X s - X 5 in n 
V-0 ■ no = on Oil, 



Xx. 



i/> = 



and choose £ := V^> in Step 3. We observe that \Xg— X§\ < 1 and \VX$\ < C(Q)5 1 . 
By standard elliptic estimates we conclude that 



c 2 (n) < -^C(Q). 



Moreover, we compute that 



\X-X S \W(a) < C(Q)5\1 + J d\VX\ 



and 



Xx 



1 

W\ 



x, < 



\n\ 



Xx 



m 



(3.27) 



(3.28) 



(3.29) 
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Therefore we deduce the estimate 



/ Xdiv£dx = [ X(X 5 -X 5 )dx = (1 - X s )m + [ (X S -X)X 
Jn Jn Jn 

m -C{Q)6(l + J d\VX\^j 



> 1 



\n 

> c(m Q ,tt), 

for 5 = 5 (m Q ,Q)(l + j Q d\VX\) . Further we compute that 



(3.30) 



n 
< 

< 



NX 



\vx\ 



d\VX\ - [ Xdiv{& )dx 
Jn 



\c 2 (n) 

5 



[ d|V^|+2||/x ||Hi ( 
Jn 



(n)\\W\\c 2 (n) 



d\VX\+C(n,n)-\\Vfi \\ L , in) 



< C{m ,n,n)(l + J d\VX\)(J d\VX\ + ||V//, 



o\\L 2 (n) 



(3.31; 



where in the last two steps we have used (3.27) and Poincare's inequality. We now 
obtain (3.14) from (3.26), (3.30), and (3.31). The estimate (3.15) follows again from 
Poincare's inequality. □ 



Proof of Proposition 3.1 We construct iteratively time-discrete solutions v h ,X h . 
First set 



v\t) := v , X\t) := X 



for -h<t<0. 



Given functions v h (t- h), X h (t- h) for t G (kh, (k + l)h] C [0,T], k G N , Lemma [3^3 
yields a solution X h (t) G W (mo) (0; {0, 1}), ^(t) G H^Q), \ h (t) G R that satisfy 
fl3~3j ), Q and 

I d\VX h (t)\ + ^\\V^(t)\\l 2{n) < I d\VX h (t-h)\ + hv h (t-h)\\l 2{n) , (3.32) 
\\^(t) + X h (t)\\ mm 

< C(mo,n,fi)(l + J d\VX h {t)\^J d\VX h {t)\ + || V/(t)|| L2(Q) ) . (3.33) 

the existence of v h {t) G Hl a {yt) that satisfies (3.1) 



Then we deduce from Lemma 
and the estimate 

1.. K 



3.2 



v h (t)\\l 2{n) + h I u(X)\Dv\t)f 
Jn 

< \\\v\t - h)f» {a) + fr||V/^)IU 2 (n)||A*)iU^) 

< \\\v h (t - h)\\l 2{n) + ^\\V^(t)\\l Hn) + h\\v h (t)\\l 2{n) . 



(3.34) 
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By construction v h , X h are constant on each subinterval (kh, (k + l)h] C [0, T], k G M. 
Summing (3.32), (3.34) for t^ = kh, k = 1, [t/h] we obtain that 



I r rh[t/h] p | 

^\\v h (t)\\l 2{n) + / d\VX h (t)\ + / / (u mm \Dv h (r)\ 2 + -\V^ 

<^IWI| a( n)+ / rf|V^ | 



"Ml 2 



da; dr 



+ 2 



h[t/h] 

h Jn 



v h (T)\ 2 dxdr 



1 fh[t/h] 



4 



Vu£(t)| 2 dx dr. 



(3.35) 



Using Gronwall's Lemma we deduce (3.5). Finally, (3.6) follows from (3.5) and 
fl3~33j). □ 



4. Passing time-discrete approximations to a limit 

We first show strong compactness of v h , X h . To this end, we will apply the following 
theorem by Simon [181 Theorem 6] 

Theorem 4.1. Let X C B C Y be Banach spaces with compact embedding X ■=— > B. 
Let T > and let J 7 be a bounded subset of L q (0,T; X) , 1 < q < oo. Assume that 
for every < ti < t 2 < T 

sup ||r s / - /|Ui( tl , t2; r) ^ as s -> 0, (4.1) 

where r s f(t) := /(i + s) /or every i G (0, T — s). Tnen JF zs relatively compact in 
L p (0, T; 5) /or every p G [!,<?). 



First of all, because of (3.1), r /l (. — h) ■ Vv h = diy(v h (. — h) <g> v h ), and since 
v h (. -h)®v h e Lt(0,T;L 2 (fi)) is bounded, d~ h v h G Ll(0,T; H~\tt)), h G (0,1) is 
bounded. This implies that 

\\r kh v h - « h || LS < V^z W-) < Ofc/i (4-2) 



for all fc G N such that fc/t < T. Therefore (4.1) holds for T = {v h : h G (0, 1)} and 
y = H~ x since r fe is p iecewise constant. Moreover, since v h G L 4 (0, T; id2 (f2)) is 
bounded, Th eore m 4.1 implies that (v h ) 0<h< i is relatively compact in L 2 (Q T ). 

Similarly, Q implies that d~ h X h G L 2 (0,T; H~\Q)), h G (0,1) is bounded. 
Moreover, since X h G L°°(Q,T; BV(tt)) is bounded and BV(fi) ^ L x (fi) compactly, 
we obtain that X h G L 1 (fir), < h < 1, is relatively compact. Since || A^H^oo^) = 
1, (X h ) 0<h<1 , is relatively compact in L P (Q T ) for every 1 < p < oo. 

As a corollary we obtain the compactness of time-discrete approximations. 

Proposition 4.2. Choose a sequence h — > and Zet (v h , X h , fi h , X h ) denote the time- 
discrete solutions constructed in Proposition 3.1 Then there exists a subsequence 
h k -> (jfe -> oo) and n G L 2 (0, T; F 1 ^) fl L°°(0, T; L 2 (fi)), /i G L 2 (0, T; rl 1 ^)), 
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X G L°°(0,T;BV imo) (Q;{0 : 1})), and A G L 2 (0, T) such that 



v k - 




->oo 


V 


in 


v k - 




->oo 


V 


in 


Xk - 


-*k- 


-+oo 


X 


in 


^k ~ 




->oo 


/i 


in 


Afc, - 




->oo 


A 


in 



for all 1 < p < oo, 



(4.3) 
(4.4) 
(4.5) 
(4.6) 
(4.7) 



where (v k , X k , fi k , X k ) := (v h " , X h " , fi h " , X h ") , G N. 



4.1. Convergence in (3.1) and (3.3). We first verify the equations in the bulk. 
Proposition 4.3. Letv,p,X be the limits obtained in Proposition ^ ■ 2 Then 

I (-vd t <p + v- Vvp + u(X)Dv : Dip) d(x,t) 
Jn T 

XV p ■ (p d(x, t), 



<p(0, x) ■ v Q (x) dx = - / XV p ■ (pd(x,t), (4.8) 

Jn T 

holds for all up G L7°°([0, T}; Cg^(fi)) with ip\ t=T = and 

/ V/i-V(d(x,t) = / a t C^ + div(C'u)A'd(x,t) + ({0,x)X {x) dx (4.9) 
/ioWs /or all ( G C°°([0,T] x H) wit/i Clt=T = 0. 

Proof. If we test in Q with </?(•, t), where p G (7°°([0, T]; Cg^(fi)) with p\ t=T = 0, 
and integrate over t G (0,T), we deduce that 

J ( - ^<9+^ + v h {. - /i) ■ Vv h <p + u(X h )Dv h : Dp) d{x, t) 

-\ f I <p(x,t)v (x)dxdt= - ! X h Vp h ■ <p d(x, t) 
h Jo Jn Jn T 



for all sufficiently small h > 0, where we set ip(t) = for t > T. By (4.3), (4.4), (4.6) 
we can pass to the limit h — > in this equality and obtain 

J ( - v h d t <p + v h ■ Vv h + v(X)Dv : Dp d(x, t)(pj d(x, t) 

ip(0,x)v o (x)dx = - XV p- p>d{x,t). (4.10) 



Similarly we obtain from (Q that for all C G L7°°([0,T] x SI) with (\ t=T = 
V/j h -V(d(x,t) 



n T 



f {x h d+ h C + (v h x h )(. - h) ■ VC) d(x,t) + [° [ X 

Jn T J-h Jn 



'o(x)((x, t) dx dt 



holds and again we can pass to the limit in this equality and obtain (4.9). 



□ 
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4.2. Convergence in the Gibbs— Thomson law. The main difficulty in passing 
the approximate solutions to a limit is the convergence in the Gibbs-Thomson con- 
dition ( |3.4 ). In particular, we cannot exclude that parts of the phase boundary 
d*{X h (-,t) = 1} cancel in the limit h -> 0. To overcome such difficulties we consider 
the limit of the phase boundaries in the sense of measures and use varifold theory. 
For the definition of varifolds and weak mean curvature for varifolds we refer to [19J. 
Let ^ := d\VX h (-,t)\ denote the surface measure of the phase interface 

d*{X h (;t) = l}, 

$ h M := f V d\VX h (-,t)\ for rj G C ({1) (4.11) 
Jn 

and let n h (t) denote the inner normal of d*{X h (-,t) = 1}, i.e., 

h , . VX h (-,t) . . 



which is well-defined for ft^-almost all x G d*{X h (-,t) = 1}. By fl3~12| the first 
variation of d 1 } is given as 

80* (ri) = [ (divr]-n h {-,t)-Dr}n h {;t))d\VX t h \ 
Jn 

X h (t) div (^(;t) + X h (t))rj)dx 



for all 7] G C£(Sl,R d ). 

We will prove that for almost all t G (0,T) the phase boundary d*{X(-,t) = 1} 
has a generalized mean curvature in the following sense. 

Definition 4.4. Let E C VL and X E G BV(O). If there exists an integral (d — 1)- 
varifold $ on Q such that 

d*E C supp(tf), 
has weak mean curvature vector H$, 
H# G Lf oc (tf), s>d-l,s>2 

then we call 

H := H$\q* e 
the generalized mean curvature vector of d*E. 

This definition was justified in [16j, where it is shown that under the above con- 
ditions if is a property of E and independent of the choice of Moreover, for 
any C 2 -hypersurface M C M. d the mean curvature H of d*E coincides Tt d ~ l -almost 
everywhere on M H d*E with the mean curvature of M. 

Lemma 4.5. Let s = 4 if d = 3 and 1 < s < oo arbitrary if d = 2. Then for almost 
all t G (0, T) the phase boundary d*{X(-,t) = 1} has a generalized mean curvature 
H{t) G L s (d\VX h (-,t),M. d )\, 

\H(-,t)\ s d\VX(;t)\ < Climinf \\/j, h (;t)\\ m{n) . (4.12) 
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Further, H(-,t) determines the limit of the first variations 8^: For any subsequence 
hi — ► (i —>■ oo) of h — > such that 

limsupll//^-^)!!^) < oo (4.13) 



and for all rj G Cq(Q, M. d ) we obtain 



60* {r,) 



-H(t)-£d\VX(;t)\ 



as i — > oo. 



(4.14) 



Proof. By Fatou's Lemma and (3.5) we deduce that t i— > liminf^o H/^Oj Oll^fn) 
belongs to L 2 (0,T) and that the right-hand side of (4.12) is finite for almost all 
t G (0,T). In the following let t G (0, T) be such that liminf/ l _ s . H/^-j t)||fri(n) is 
finite. 

Since J n d\VX h (-,t)\ is uniformly bounded by ( |3.5 ) and recalling (4.13) we can 
extract a subsequence (not relabeled) hi — > 0(i — > oo) such that 

— >■ $ t as Radon measures, (4-15) 

fi hi (-,t) w t in if 1 ^), (4.16) 

for a Radon measure ?9 4 on Q and w 4 G H 1 ^). We then deduce from [TTJ Theorem 1.1] 
that 

$ t is an integral varifold, 
— > as varifolds, for a subsequence /ij — > (i — ► oo), 
fit has weak mean curvature H$ t G L s ($t), 
|VAT(-,t)| < t , 



and that 



holds i?t-almost everywhere, with 



Ha 



n{-,t) 



\VX(;t)\ 





w t n(t) 

on d*{X(-,t) = 1}, 
elsewhere. 



According to [16] H(-,t) := H^ t \ s*{x( -,t) =\} is a property of and independent 

of the choice of subsequence in (4.15), (4.16). Moreover, due to [T7J Theorem 1.2], 
we have 

H(-,t) = -^-almost everywhere in {9 d ~ 1 ($ t , •) ^ 1}, 

where 9 d ~ 1 ($ t ,-) denotes the (d — l)-dimensional density of $ t , cf. [19]. Since the 
first variation is continuous under varifold-convergence, we then obtain that 

lim Stipto) = 6MV) = [ -H^-rjd&t = [ -H(-,t) ■ rj d\VX(-,t)\. 
i-* 00 Jn Jn 

□ 

We still have to relate the generalized mean curvature H(-,t) that we obtained for 
almost all t G (0, T) with the weak limit \i of // in L 2 (0, T; H\n)). 
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Lemma 4.6. For all £ G L 2 (0, T; Cg(fi; R d )) 



[ -H(-,t)-£{-,t)d\VX{-,t)\dt= I X(x,t) div (fi(x,t)£(x,i)) d{x,t). (4.17) 
Jo Jn T 



In particular, for almost all t G (0,T) ; 

H(;t) = fj,{;t)n{;t) 

holds 7t d ~ l -almost everywhere on d*{X(-,t) = 1}. 

Proof. Define for all t G (0,T) such that H{-,t) G L s {\VX(-,t)\) exists 



(T(t),^) := / -H(x,t)-^d\VX(x,t)\ 
Jn 

for all V e Cft^R*). Then 

|<r(t),V>| < C|^||co(a)||^(-,*)llM( d |v^|) 



and we deduce from gj§ that T G L 2 (0,T; C (fi)*). 
Similarly we define for ^ G Cq(^) 

(!*(*), V>) := Wf(V)- 



From (3.4) we deduce that 

< C||Ai fc (*)ll^(n)lie(*)llcS(n), 



and 



l^ ft |U 2 (0,T;C3(n)*) < C||^ /l |U 2 (0,T;H 1 (n))||^||L2( 0>T;C -i(Q)). 



Moreover, by (4.6), (4.7) there exists a subsequence h — ► such that 



lim / (T ft (t),e(-,t))rft 



o Jn 



# V ■ (£//) da; dt. 



For a > we now define functions T% : (0, T) -> Cg(n; M d )*: 

J<7?,V>> if ||Ai*(-,t)||fli ( iD < 



CC(t),V>> 



a. 



\(T(t),V> else. 



(4.18) 



(4.19) 



(4.20) 



Fix an arbitrary £ G L 2 (0, T; Cq (fl; R ')). Then we deduce from Lemma 
almost all t G (0, T) 

(T£(t), £(-,*)} - [ -H.a-,t)d\VX(.,t)\. 



4.5 



We also see from (4.21) that 



|(T*(t),£(-,t)>| < ||e(-,t)||cS(n^)C(«+l|r(t)||cb(o) 



(4.21) 
(4.22) 

(4.23) 

(4.24) 
that for 

(4.25) 
(4.26) 
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which gives by (4.20) an L 1 (0, T)-dominator for the left-hand side. By (4.25), (4.26) 
and Lebesgues Dominated Convergence Theorem we deduce that for all a > 



(T(t)£(-,t))dt as/i->0. 



(4.27) 



Next, consider the sets A := {t G (0, T) : t) Wmtn) > a } an d observe 



T 



< 



(T h (t)-T^t),a-,t))dt 

ll^('' i )llc , i(Q;IR d )^ 



< 



rph I 



<7*(f)-T(f), £(.,*)) 



dt 



L 2 (0,T;C*l(O;] 



)*) + ll^llL 2 (0,T;C*o(n;l 



By (|420[ and gp) \\T h \\ L 2 (0 ^ c i 
formly in h > 0. Since 



(o,T;C5(n^«»)*) and II t IIl 2 (o,t ; c ( 5(W*) are bounded uni- 



l^ ft | — a 2 ll/ /A HL 2 (0,T;.ff 1 (n)) 



< 



we end up with 



/ (T h {t)-T%(t),t{;t))dt -»• asa^oo 



uniformly in ft, > 0. Thus we obtain from (4.25), (4.28) that 



X div(^) tZ(a;, t) = lim / A? div(//£) t) 



(4.28) 



(4.29) 



lim 



T 



(£(-,t),T{t))dt, (4.30) 



which proves (4.17). Since no time derivative is involved, we deduce that for almost 
all t G (0, T) and all £ G C^fi; M d ), 



A-(-,t)div (/i(-,t)e)^ = (T(t),e) 
holds. The Gauss-Green theorem [7J Theorem 5.8.1] and ( |4.19[ ) yield 



/ H(.,t).£d\VX(.,t)\, 



(i(.,t)n(.,t)-Zd\VX(;t)\ 
with n{-,t) = ™ x \% on d*{X(-,t) = 1}. This finally proves pl8| |. 

Appendix A. Sharp Interface Limit 



□ 



Here we discuss the relation between (!.!)-( 1.10 ) and its diffuse interface analog 



(1.11 )-( 1.17). First we consider the corresponding energy identities. For the Navier- 
Stokes/Mullins-Sekerka system we recall that by (1.18) every sufficiently smooth 
solution of ( 1.1 )-( 1.10) satisfies 

d 1 



dt2 



[ \v{t)\ 2 dx + K^-H d -\T{t)) = - [ v{X)\Dv\ 2 dx - m [ \V^\ 2 dx. (A.l) 
Jq dt Jq Jq 
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On the other hand, every sufficiently smooth solution of (1.11)- (1.17) satisfies 
d 1 



l v(t)\ 2 dx + ^-E £ (c(t)) = - [ u(c)\Dv\ 2 dx - m [ |V/i| 2 cfe, (A.2) 
dt 2 J n dt J n J n 



where 



EJc) 



[ |Vc| 2 rfa; + - / f{c)dx. 
in £ Jn 



Moreover, by Modica and Mortola [14J or Modica [13J, we have 

E £ — > £ ^ V w.r.t. L 1 -r-convergence, 

where 

tzH d ~ l (d*E) if u = 2Xe — 1 and E has finite perimeter, 



V{u) 



-oo 



else. 



Here k = /„ v / 2f(s)ds, /: R — > [0, oo) is a suitable function such that f(s) = if 
and only if s = 0, 1, and d*E denotes the reduced boundary. Note that d*E = dE 
if E is a sufficiently regular domain. Therefore we see that for constant m > the 



energy identity (A.l) is formally identical to the sharp interface limit of the energy 
identity (A.2) of the diffuse interface model (1.11)- (1.17). In contrast, if we would 
choose m = m £ ^ £ ^o in (1.11 )-( 1.17), we formally obtain in the sharp interface 



limit the energy identity of the classical two-phase flow (1.1 )-( 1.10) 



Now we will adapt the arguments of Chen [5] to show that as e 
m(e) — > 



and if m 



m > solutions of the diffuse interface model (1.11 )-( 1.17) converge to 



varifold solutions of the system (l.l)-(l.lO), which are defined as follows: 



Definition A.l. Let v G L^.(Q), E C £1 be a set of finite perimeter, and let 
Q = Q x (0, oo). Then (v,E, V) is called varifold solution of (1.1)-(1.10) with 
initial values (vo,E ) if the following conditions are satisfied: 

(1) v e L°°(0,oo;L2(fi)) n2/ s (0 > oo;£#(ft) d ), fi G L 2 OC ([0, oo); H\Q)) andVfi G 
L 2 (Q). 

(2) E = [J t>0 Et x {t} is a measurable subset of Q x [0, oo) such that X% G 
C([0, oo); L\®)) n L°°(0, oo; BV(Q,)), J Q X E (t)dx = m for all t > 0, and 
X E \t=o = X Eo in L x (fi). 

(3) V is a Radon measure on Q x Gd-i x (0, oo) such that V = V t dt where V 1 
is for almost all t G (0, oo) a general varifold on Q, i.e., a Radon measure on 
Q x Gd-i- Moreover, for almost all t G (0, oo) V 1 has the representation 



ip(x,p)dV t (x,p) = / b i(^( x ,Pi( x )) dxt (3c 
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for all ip G C (Q x Gd-i), where A* is a Radon measure on Q, b\,. . . ,b* N and 
p\, . . . ,p l N are measurable R- and W 1 -valued functions, respectively, such that 



< b\ < 1, V^fo* > 1, ®p\ = I /J -almost everywhe 

i=\ i=l 

< — X t -a.e. with k— a/2/(s) ds. 
K Jo 



X f 

(4) Moreover, 

/ {-vd t (p + v ■ Vvp + v{X E )Dv : Dip) d(x, t) 
Jq 

— / <p\t=o ■ v dx — — k / X E V fx • if d(x,t) (A.3) 
Jn Jq 

holds for all ip G C°°([0, oo); Cg^(fi)) uritfcsuppp C fix [0,2*1 for someT > ; 

m / Vfi-Vtjjd(x,t) = / ^(9i^ + div(^))d(x,t) + / V|t=o^E ^ ( A - 4 ) 
Jq Jq Jn 

holds for all G C°°([0, oo) x O) ratt supp^ Cflx [0,T] /or some T > 0, 
as u>eW as 



(^,div( W )) n = (oV*,7?) := / {I -p®p):Vr ] (x)dV t (x,p) 

for all f] G Cq(Q x 
(5) Finally, for almost all < s < t < oo 

^IK*)lli" ( n) + A'(n) 

+ T f (v(X ET )\Dv\ 2 + m\V^\ 2 ) dxdr<hv(s)\\l 2{u) + X s (U). (A.5) 

Here and in the following / shall satisfy / G C 3 (1R), f(s) > and f(s) = if and 
only if s = 0, 1 as well as f"(s) > c (l + |s|) p ~ 2 if s > 1 — c and if s < c for some 
c > 0,p > 2. Then F(s) := f(^) we will satisfy the assumption in [3]. We will even 
assume that p > 3, which we will need in the following to estimate v e - Vc e = div(t> e c e ) 
uniformly in L 2 (0,T; H~ 1 (Q). One can choose e.g. f(s) = s 2 (l — s) 2 . 

For the following we denote 

e £ {c) = S -\Vc\ 2 + e- 1 f{c). 

Theorem A. 2. Let Vt C W 1 , d = 2,3, be a smooth bounded domain and let v G C°(R) 
with v[s) > Vq > for all s£i Moreover, let initial data (t>o,e> c o,e) G L 2 (Q)xH 1 (fl) 
with p| J Q co, £ dx = c G (—1, 1) 6e azwen iaai satisfy 

\ f \v , £ \ 2 dx + E e (c , £ ) <R (A.6) 
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uniformly in e G (0,1] for some R > 0. Finally, let (m e ) E6 (o,i] with m £ 



m > 0. Consider now (l.ll)-(l.ll) with m replaced by m £ and let for every e G 



(0,1] (v e ,c e ,fi E ) be weak solutions in the sense of [TJ Definition 1]. Then there is a 
subsequence (ek)keN such that \ as k — ► oo such that: 

(1) There are measurable sets E C fl x [0, oo) and EoCfi swc/i i/ioi c £fc — >k->oo Xe 
almost everywhere infix [0, oo) and in C«([0, T]; L 2 (f2)) /or allT > as well 
as co tEk —>k~^oo —1+2Xe almost everywhere in fl and Xe = XeU=o in L 2 {fl). 

(2) There are \i G L 2 OC {[0, oo); H\fl)),v G L 2 (0, oo; H^(fl)), v G such that 

/i in L 2 (0, T; if 1 ^)) /or a// T > 0, (A.7) 
W« mL^O.txjF 1 ^)), (A.8) 
V-** »n (A. 9) 

, . . . , n on fl x [0, oo) such that 



(3) There exist Radon measures A and \j, i,j 
for every T > 0, i,j — 1, . . . , n 



e £fe (c £fc ) dx dt 
£kd Xi c £k d Xj c £k dxdt 



k— >oo 



A 



k— >oo A «7 



mM{fl x [0,T]), 
in X(flx [0,T]). 



(A.10) 
(All) 



There exist a Radon V = V* dt on fix Gj-i x (0, oo) such that (v, E, /i, V) is 
a varifold solution of (1.1)-(1.10) with initial values (v ,E ) in the sense of 



dt 



Definition \A.l\ where 

f (SV*,ri) dt = [ [ V?7 : (dXI - (dAy)& =1 ) 
Jo Jo Jn 

for all 7]e Cg(fi x [0,T];M d ). 
First of all, by the definition of weak solutions in [I], we have 

-v e d t if> + v £ ■ Vv £ 4> + v{c £ )Dv £ : DijS) d(x, t) 



(A.12) 



v ip\ t=0 dx 



c e V/i e • if) d(x, t) 



(A.13) 



for all i) G L7°°([0, oo) x fl) d with div^ = and ip(t) = for t > T for some T > 0, 
as well as 

m £ I V// E • Vipd(x,t) = / c e (d t ip + div (ipv e )) d(x,t) + / c 0i£ (p\ t=0 dx (AAA) 
Jq Jq Jn 

/ fi e ipd(x,t) = / f(c e )<pd(x,t)+ / Vc e • V<pd(x,t) (A.15) 

for all ip G C°°([0, oo) x fl) with supp ^Cfix[0,T] for some T > 0. Moreover, from 
the energy inequality in [TJ Definition 1] we obtain 

l\\ve(t)\\ 2 L 2 {n) + E £ (c £ (t)) + [ [ (u(c £ )\Dv £ \ 2 + m £ \Vfi £ \ 2 ) dxdt<R 
* Jo Jn 
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for all t G [0, oo). Therefore we have 

IM-,f)||zri(n) < C (E £ (c £ (t)) + ||V y u £ (.,t)|| L2( n ) ) (A.16) 

for all t > and < e < Eo for some C, Eq > due to [SJ Lemma 3.4]. 

Hence there exists a subsequence Ek \ as k — > oo such that ( A.7)-( A~9~| ) holds. 
Moreover, using (1.11) and the Lemma by Aubin-Lions, one easily derives that 
v £k — > fc ^oo v strongly in L 2 (Q x (0,T)) for all T > and v £k (t) —^k^ v(t) strongly 
in L 2 (Q) for almost every t > 0. 

Using the assumptions on / we further deduce that 

/ \c £ (t)\ p dx < C{1 + R), [ dist(c, {0, 1}) 2 dx < CeR 
Jn Jn 

uniformly in e G (0, 1], for almost every < t < oo and t = 0. Now we define as in 





w £ (x,t) = W(c £ (x,t)) where W(c) = / v/2/(s) (is, /(s) = min(/(s), 1 + \sf). 

Jo 

Then (Vio e ) £e (o,i] is uniformly bounded in L°°(0, oo; L 1 (J1)) since 

/ \Vw £ (x,t)\dx = / J2f{a)\Vce{x,t)\dx< I e £ (c £ (t)) dx < R. (A.17) 
in in Jn 

Moreover, we have for all ui,U2 G 1 

ci|«i - u 2 \ 2 < \W{ui) - W(u 2 )\ < c 2 \u! - u 2 \ 2 (l + + \u 2 \) (A.18) 
which again follows easily from the assumptions on /. Now we obtain: 
Lemma A. 3. There is some C independent of £ G (0,1] such that 

Proof. The proof is a modification of [3 Lemma 3.2]. Therefore we only give a brief 
presentation, describing the differences. 
For sufficiently small rj > let 

c v e (x,t) = p(y)c £ (x -r]y,t)dy, xefl,t>0, 

JBi 

where p is a standard mollifying kernel and c £ is extended in an ^-neighborhood of 
fl as in [5j Proof of Lemma 3.2]. Then one obtains 

\\Vc%,t)\\LHu) < Ctj-'Wc^WlHu) < CV 1 
||c?(.,t) - c £ (.,0ll^ ( a) < C77||V^ e (.,t)|U 1( n) < C77 
for all sufficiently small rj > 0, cf. P2 Proof of Lemma 3.2]. Next we use that 

(c e (;t)-c e (;T),<p) (l = - / ((V// - v £ c £ ) (s) , V(p) Q ds 



for all if G C 1 (f2) because of (1.13) in its weak form. Here 

|| V/x - v £ c e || L 2( Cx ( T]t) ) < C{R) for all < r < t < 00, \t - t\ < 1 
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since v £ G L 2 (0, oo; L 6 (£l)) and c £ G L°°(0, oo; L 3 (f2)) are bounded due to d < 3 and 
p > 3. Hence 



(c^(x,t) - c^(x,r))(c e (x,t) -c e (x,r))<ix 



((V/i - ^c £ )(s), Vc^(x, t) - V4'(x, r)) n rfs 



< C{R){t-r)> sup ||Vc2(. > s)|| £a( n)<C r (T,fl)7 7 - 1 (t-r)2 

s6[0,oo 

for all < r < t < oo, |t — r| < 1. Now, choosing r\ = min(?7o, (t — t)s), we conclude 



\c £ (.,t)-c £ (.,T)\\ 2 L2{n) <C\t-r\^ 



for all It - r| < 1 



from the previous estimates. Thus c £ G C» ([0, oo); L 2 (f2)). Using (A. 18), one derives 
U7 e G CsflO.oa);!^)) as in 0. 



□ 



Note that the previous lemma and (A. 17) implies that (w £ )o <£ <i is bounded in 
L°°{0,oo;BV{Q)). 

Lemma A. 4. There is a subsequence Ek \ as k — > oo, a measurable function S(t), 
t G (0, oo) and a measurable set E C Q x (0, oo) stzc/i t/iat 

(%) E ek (c Ek {t)) -> £(t) for almost allt>0. 

(2) w £k — > almost everywhere in Q x (0, oo) and zn Cs ([0, T]; L 1 (f2)) /or a// 
T > 0. 

f5) c £fc — >■ Afg almost everywhere in Q x (0, oo) and m Cs ([0, T]; L 2 (fi)) /or a// 
T>0. 

In addition, X E G L°°(0, oo; BV(fi)) n C7j ([0, oo); L 1 ^)) and E t := {x e Q : (x,t) G 



E} satisfies \E t \ = \E 



1+5 
2 



fl| /or almost all t > and 
IVAf^Kfi) < -£(t) < -£ - 

K K 



Proof. The main difference to the proof of Lemma 3.3] is the proof of convergence 
for E £k (c £k (t)). To this end one uses that 



F £ (t):=-\\v £ (t)\\l Hn) + E £ (c £ (t)), t>0, 

is a sequence of bounded, monotone decreasing functions and v £k (t) ^k^oo v(t) for 
almost all t > in L 2 (Q). The rest of the proof is identical with the proof of [3 
Lemma 3.3]. □ 



Using the previous statements one can now easily finish the proof of Theorem A. 2 



by the arguments of [HJ Section 3.5]. In particular, (A. 3) and (A. 4) easily follow from 
(A. 13) and (A. 14). It mainly remains to show (A. 12) and (A. 5). Let A, Ajj be as 
in (A.lO)-(A.ll). To show the energy estimate (A. 5) one uses that dX = d\ l dt for 
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some Radon measures A* on Q and that for almost every < t < s < oo 
A*(H) = lim E £k (c £k (t)) 



k— >oo 

ft 



k- 



-OO 



< lim I E £k {c £k{s) ) 



I I ( u ( c e k )\Dv £k \ 2 + m £k \Vfi £k \ 2 ) dxdr 

Js JQ 



+ ^[^he k {s)\\ 2 L 2 {n) ~ ^\\v £k (t)\\ 2 L 2 (n) 

I , I 



< A*(fl)- / / (u(X E )\Dv\ 2 + m\Vfi\ 2 ) dxdr+-|K a )||i 2(n) --Ht)||i a(n) , 

J S J Q 

where we have used the weak convergence of Dv £k , V/z £fe in L 2 (Q) and the strong 
convergence of v £k (t) in L 2 (Q) for almost all t > 0. Finally, (A. 12) follows in pre- 



cisely the same way as in [51 Section 3.5], where we note that the main point in the 
argumentation is that the discrepancy measure 

ffe) :=||Vc £ | 2 -^(c £ ) 

converges to a non-positive measure. The latter fact follows from [5j Section 3.5] 
since we have the same bounds on // e in 
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